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This overview provides an interface between lines of thought on rh mitigation in the 
fields of radio astronomy and signal processing. The goal is to explore the commonality of 
different approaches to help researchers in both fields interpret each others' concepts and 
jargon. The paper elaborates on the astronomers' concept of gain closure relations and 
how they may be used in a "self-calibrating" system of rfi cancellation. Further discussion 
of the eigen decomposition method in terms of rfi power and antenna gains introduces 
adaptive nulling and rfi cancellation through matrix partitioning. While multipath 
scattering appears at first glance to be fatal to methods of rfi cancellation, its effect is 
easily incorporated in frequency dependent gain coefficients under many circumstances. 
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1. Introduction 

Astronomers are designing the next generation ra- 
dio telescopes with the capabilities to solve many cur- 
rent puzzles about the origin of the planets, stars, 
galaxies and Universe itself. Since the telescope de- 
signs are one to two orders of magnitude more pow- 
erful than the best present telescopes, astronomers 
know that the new instruments are capable of mak- 
ing discoveries that will open whole new fields of re- 
search. New telescopes will not just solve old prob- 
lems, but will certainly do unexpected things that 
have been impossible in the past. 

The same advances in technology that make new 
radio astronomy possible are closely linked to the ad- 
vances in telecommunications that are beginning to 
saturate the radio spectrum. Any new radio tele- 
scope will have additional layers of signal processing 
built in from the design stages that are intended to 
enable it to see through the background radio inter- 
ference as well as possible. In this regard, modern 
technical solutions offer great hope, but in reality 
they are complex and expensive, and the resources 
(in terms of time and money) will probably not exist 
to allow radio astronomers to apply these methods 
effectively enough to keep up with the rising earth- 
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generated radio background unless there is active co- 
ordination with the telecommunications industry. 

Communications signals are fully polarized, coher- 
ent emission from carefully modulated currents in 
the broadcasting antennas; a communications signal 
is broadcast to provide maximum efficiency in use 
of a limited portion of the radio spectrum defined 
by its allocated bandwidth BW. The astronomical 
signals originate mainly in the incoherent emissions 
from huge clouds of faintly emitting atoms, molecules 
or energetic electrons whose "astronomical intensity" 
is a collective effect arising from their large num- 
bers. The radiated powers are added incoherently, 
except in a few cases, such as the coherent addition 
of radiated electric field in the ionospheres of pulsars, 
Jupiter's ionosphere, the active regions of the Sun or 
celestial masers. Astronomical signals in general are 
well represented by stochastic or gaussian noise. 

Telecommunication reception requires decoding 
the fluctuations in the electric field of the signal, so 
that the field at each time step T (where T « 1 / BW) 
reliably carries information. The information that 
astronomers extract from their signals has a more 
statistical character, such as total intensity of the 
power per unit bandwidth (called flux density) of an 
astronomical source, along with other more detailed 
statistical measures such as the net polarization of 
the signal power and the variation of the flux density 
as functions of frequency and time. For the large va- 
riety of celestial sources that fall in the class called 
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'continuum sources,' which radiate a smooth spec- 
trum over a wide range of frequencies, more power 
can be received from the source by simply increasing 
the bandwidth of the observation. 

Astronomers have not generally been interested in 
measuring the electric field strength as a function of 
time, although this is likely to change for some astro- 
nomical applications, such as rapidly time- variable 
phenomena (for example: pulsars, transient bursts, 
signals from energetic particles impinging on the up- 
per atmosphere Huege & Falcke 2003 [@]; Suprun 
et al 2003 [@]), especially with the more sensitive 
telescopes being planned for the future. The need 
for RFI mitigation has also increased astronomers' 
awareness of the need for processing their E- field sig- 
nals in real-time at full bandwidth, which is much 
more demanding and expensive than their historical 
procedures (Barnbaum & Bradley 1998 [@]). 

While there do exist frequency bands that have 
been reserved for astronomy, receivers have become 
so sensitive that low level spurious harmonics and 
weak leakage from communication and navigation 
services are strong enough to corrupt and confuse 
the astronomical signals. Furthermore, pursuing an- 
swers to scientific questions drives astronomers to ob- 
serve outside of the reserved bands. For example, 
neutral atomic hydrogen has a well known radiative 
transition at rest frequency f a = 1420.4057 MHz, for 
which there is a reserved astronomy band. How- 
ever, astronomers now wish to trace the evolution 
of hydrogen in galaxies over cosmic time, and since 
observing galaxies at redshift z « 1 would corre- 
spond to looking back in time to when the uni- 
verse was less than half its current age, radio as- 
tronomers pursue the hydrogen line to frequencies of 
/ = /o/(l + -z)~ 700 MHz or lower, where the radio 
spectrum is already in heavy use. 

Astronomers have had great success in achiev- 
ing high angular resolution and high dynamic range 
imaging. The question naturally arises, "why not 
simply image the rfi sources along with the celes- 
tial ones?" Under some circumstances, this is a suc- 
cessful approach (e.g. Perley et al 2004 [@]); how- 
ever, the character of many RFI signals - variability, 
non-gaussianity - may not be compatible with the 
necessary conditions for obtaining sufficiently high- 
dynamic range aperture synthesis. 

As in all of radio science, improvements in instru- 
mentation and signal processing are providing the 
tools for performing astronomical observations in ra- 
dio rich environments. These include: 



1. Robust receivers Modern receivers benefit 
from linearity over a large dynamic range to avoid 
generation of spurious signals through intcrmodula- 
tion within the receiving system. Since astronomical 
signals are noise-like, astronomers have historically 
processed their signals with custom built correlators 
that operate at 1 or 2 bit precision, at the cost of 
adding a modest percentage of quantization noise. 
Newer telescopes that will operate in populated ar- 
eas are designed to process data with 10 to 14 bits of 
precision, a cost only now affordable due to advances 
in low-cost digital processing electronics. 

2. Sophisticated signal recognition and 
blanking mechanisms. Historically, this was per- 
formed in the post-detection domain, where it was 
labor intensive and called "data editing." If done 
early in the receiving system chain, blanking can ease 
the large dynamic range requirements of the receiver 
backend (reducing the required bits of precision for 
example). It is always a concern that a process so 
violently non-linear as blanking may preclude the ap- 
plication of some classes of rfi cancellation for which 
linear systems are a prerequisite. These topics in 
rfi mitigationnhave been reviewed by Fridman and 
Baan (Fridman & Baan 2001 [@]) and are addressed 
by several papers at this conference. 

3. RFI cancellation A number of recent papers 
have begun to explore methods of rfi cancellation and 
subtraction(Barnbaum & Bradley 1998 [@]; Elling- 
son et al 2001 [©]; Leshem et al 2000 [@]; Briggs et 
al 2000 [@] ) . The series of papers by Leshem, van der 
Veen, Boonstra and colleagues have been especially 
important in applying the techniques of signal pro- 
cessing in the radio astronomical context (see Boon- 
stra et al 2004 [@] in these proceedings for complete 
references) . 

The bulk of this paper forms a mini-tutorial cover- 
ing several issues in rfi cancellation that were sources 
of confusion to the authors and their colleagues. The 
goal is to ease the transfer of technology between the 
astronomical and signal processing communities. 

Of course, the best starting place for radio as- 
tronomy continues to be siting telescopes in radio 
quiet environments, and astronomy will always rely 
on good will and cooperation with radio-active spec- 
trum users for communication and navigation ser- 
vices. 
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Figure 1. Astronomical telescope arrays with (1) conventional dish antennas pointed toward ce- 
lestial radio sources and (2) reference antennas to provide increased gain in the direction of known 
interfering signals. 



2. Radio Astronomical Arrays 

Modern radio telescopes operate in arrays in order 
to economically obtain sensitivity and angular reso- 
lution through the combination of the signal from the 
elements of the array (Figure 1). The radio astro- 
nomical convention is to denote antenna separations 
by a vector quantity Bij , indicating the relative posi- 
tions of the antenna phase centers for antennas i and 
j. These "interferometer baselines" are conveniently 
measured in wavelengths B\ for which the effective 
angular resolution of each interferometer pair in the 
array is 6ij = 1/B\^j radians. Definition of a carte- 
sian coordinate system (u-v-w) with its w axis point- 
ing in the direction of a radio source allows B\ to be 
specified by its components u — B\, u , v = B\. v , and 
w = B\ w . An array of elements at fixed locations 
will experience systematic variation of u, v and w as 
the Earth rotates. The w component specifies the 
relative propagation delay between two elements of 
the array tij — Wij\/c, while the u and v components 
appear in the key equation of aperture synthesis, 



V{u,v) = — p B {x,y)T b {x,y)e a ^ x +^dxdy (1) 



relating the "fringe visibility" V(uij,Vij) sensed by 
each interferometer pair i,j in response to the an- 
gular brightness temperature distribution of the sky 
Tb(x,y) modulated by pB(x,y), the primary beam 
gain function of the array elements. The strength 
and phase of the complex visibility result from the 
fourier integral of pB(x,y)Tb(x,y) over the celestial 
source. Here, k is the Boltzmann constant, and x 
and y are angular coordinates measured in the plane 
of the sky. While this integral and its inversion form 
the essence of the aperture synthesis technique, a full 



description of the method requires entire books (e.g. 
Thompson et al 1986 [@]; Taylor et al 1999 [@]). 

Interferometer arrays measure a V m (u, v) for each 
baseline through the cross correlation dj of the sig- 
nals Si and Sj from two elements: 

V m (uij,Vij) = <g l s l g*s*> (2) 
= 

The signals entering the correlator have been scaled 
by the complex gain factors gi that quantify the 
differences in the signal paths from one array ele- 
ment to the next, and * indicates complex conju- 
gation. In general, V m (u,v,f), gt(f) and Sj(/) are 
functions of frequency, and Cij(f) is the observed 
cross power spectrum. Provided the gains are con- 
stant over the integration time of the cross corre- 
lation product in Eqn. 2, the expression becomes 
Vij = V m (uij,Vij) = g l g* < s,s* >. 

3. Gain closure 

The concept of "closure relations" becomes clear 
when combinations of observed correlation products 
are formed: 

q _ VjjVkn _ 9i9j < f£f| > 9k9n < S k S* n > 

V in V k j gigl < s l s^ > g k g* < s k s* > 

CijCkn r n \ 
= TT-^- (3) 

The quantity Q is a ratio of true visibilities for the ra- 
dio source, regardless of whether the individual com- 
plex gains gi are known accurately or not. 

The special case of an unresolved point source is 
of interest, since substituting a <5-function for TJ,(0, 0) 
in Eqn. 1 yields a visibility function that is constant 
for all u, v, implying that Q = 1 and the < s^s* > 
are equal for all so that the right hand side of 
the equation reduces to an identity relation among 
the gain coefficients. 
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The gain closure relations form the motivation for 
approaches to rfi mitigation since rfi sources can be 
treated as point sources and the cross correlations 
between rfi reference antennas and the astronomical 
signals provide the necessary information for com- 
puting the rfi contamination without precise knowl- 
edge of the gains of the element sidelobes through 
which the rfi enters the telescopes. These methods 
are in essence "self-calibrating." 

4. Cancellation 

Most rfi cancellation schemes involve first identify- 
ing the interferer, followed by deducing the strength 
of its contamination within the astronomical data 
stream, and finally, precise subtraction of interferer's 
contribution. 

Fig. 2 illustrates one example of how that might be 
accomplished using the signal / from an rfi reference 
antenna that is cross correlated with the astronom- 
ical signal Sa in order to deduce the inter-coupling 
of the two signal paths. Several variations of imple- 
menting this scheme through the gain closure rela- 
tions of cross power spectra were explored by Briggs 
et al 2000 [@]. Other schemes are implemented in 
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Figure 2. Example of interference cancellation scheme. 
In this case, there are two signal paths: (1) the astronom- 
ical signal Sa entering the main antenna but corrupted 
by interference 7, and (2) a reference antenna to pick up 
the interfering signal apart from S. The goal is for the 
signal processing to obtain a suitable signal for subtrac- 
tion from the astronomical signal path to leave only an 
undistorted S. Techniques discussed here use cross cor- 
relation between Paths A and R to determine the filter 
coefficients. 



the time domain (e.g. Barnbaum & Bradley 1998 
[@]; Kesteven et al, 2004 [@]). 

Separate rfi reference antennas serve two purposes: 
(1) They improve the signal-to-noise ratio on the in- 
terfering signal, in order to minimize the addition of 
noise to the astronomical observation during the ap- 
plication of rfi subtraction, and (2) with increased 
gain toward the interferer but decreased gain toward 
the celestial source, they avoid inadvertent subtrac- 
tion of the astronomical signal along with the inter- 
ference. It is important that the interference-to-noise 
ratio be better in the reference channel than the as- 
tronomy channels (Briggs et al 2000 [@]), but this 
is generally true since rfi enters the astronomy tele- 
scopes through weak sidelobes (roughly at the level 
of an isotropic antenna ^0 dBi), while simple refer- 
ence antennas can deliver 10 dB or more gain). 

The gain information contained in the cross power 
spectrum Car{J) between (qaSa + 9 Ail) and g RI I 
allows the rfi contamination of the voltage spectrum 
9AiI(f) hi the astronomical data stream to be esti- 
mated from the reference voltage spectrum gnil: 



9aiI = 



Car _ 9ai9 R iII* 



ah 1 * 



ah 1 * 



(4) 



Astronomers commonly measure the power spec- 
trum of the signal from the telescope Pa{J), which 
will have both power from the celestial source and rfi 
components Pa(J) — \9aSa\ 2 + \dAiI\ 2 - In this case, 
an appropriate correction is preferably derived from 
two reference antennas labeled Rl and R2: 



\9AiI\ 2 = 



9Al3h\I\ 2 d*AldR 2 \ 

9*ri9 R2 \ i \ 2 



(5) 



The preference for dual reference signals arises from 
the nature of the noise in the cross power spec- 
trum in the denominator of Eqn. 5; the noise term 
<NfnN R2 >(f) is a complex cross power spectrum 
that averages toward zero as integration time in- 
creases, unlike the noise <\Njn\ 2 >(f), which is real 
and positive and which would bias the result of 
Eqn. 5 if it were based on a single reference sig- 
nal. When the subtraction is performed after the 
calculation of power spectra in this way, it is often 
described as "post-correlation" cancellation. 

Similar estimation for the contamination of the 
V m (u, v) measurement of an interferometer pair i,j 
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of an array (Fig. 1) is 
gAii9Aji\I\ 2 = 



C AiR\C* AiR1 

CI 



(6) 
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with dual reference antennae. 
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Figure 3. Application of post-correlation gain-closure 
cancellation. Observations in the 50cm digital TV band 
from the Parkes Telescope using the CPSR2 recording 
system. The 64 MHz wide bands were analyzed in 16384 
channel spectra. Top panels show reference signals from 
two yagi reference antennas aimed at the broadcasting 
tower. The astronomy total power spectra for the two 
polarizations of the Parkes Telescope are shown as 1) un- 
corrected (central panels) and 2) corrected (lower pan- 
nels). 



Figure 3 gives an example of applying the rela- 
tions in Eqn. 5. These Parkes observations covered 
a 64 MHz band centered on 676 MHz with 16384 
spectral channels (3.9 kHz resolution). The goal of 
the program is the observation of neutral hydrogen 
clouds at reshift z ~ 1 in absorption against high z 
radio galaxies. Two yagi reference antennas were ori- 
ented to point at the broadcasting tower for digital 
TV station channel 46 (652-659 MHz). The cancella- 
tion algorithm was applied to ~0.1 sec integrations, 
followed by averaging for the full 16 sec scan. Af- 
ter cancellation of the digital TV, the noise level and 
weaker rfi contamination remains at the same level as 
the surrounding frequency bands; only those interfer- 
ers in the yagi reception pattern have been canceled. 
(Kesteven et al, 2004 [@] give more details on the na- 
ture of the digital TV interference in their discussion 
of a real-time adaptive canceler.) 

Figure 4 provides a closer look at the band con- 
taining the digital TV as a function of time through 
the 16 second scan - both before and after cancel- 
lation of the rfi signal. The noise level should be 
compared to that in a relatively "clean" band away 
from the TV (in the top panel of the figure). The 
clean band also contains an image of the peculiar 
drifting interferer, one of whose harmonics appears 
to the left of the TV band in panels B and C. The 
weak rfi spikes that remain throughout the spectrum 
(1) are narrower than the expected HI signals, (2) 
appear as emission features rather than absorption 
lines, and (3) take up a minor fraction of the spec- 
trum as evidenced by the broad expanses between 
the narrow lines. This is encouraging for the success 
of redshifted HI studies. 

5. Multipath delays 

At first glance, multipath propagation or scatter- 
ing of the interfering signal on the way from the 
transmitter to either the astronomical or reference 
antennas would appear to corrupt the interfering sig- 
nal - in different ways along the paths to the different 
receiving systems - rendering it impossible to extract 
an appropriate reference signal. Fortunately, multi- 
path is not a problem provided the paths are stable 
on sufficiently long time scales. 

This can be easily seen by noting that the inter- 
fering signal received by an antenna i r (t) is the sum 
of different versions of the signal i (t) broadcast by 
the transmitter at the source: 
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A. RFI Free Region 




B. Digital TV band - uncancelled 
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Figure 4. Dynamic spectra for A) a spectral region free of digital TV contamination, B) the 
spectral band containing TV channel 46, prior to cancellation, and C) the channel 46 band after 
cancellation. 



ir{t) = ^2h a (t)*i (t-T a ) 

a 

Here h a (t) is the impulse response experienced by 
the interfering signal along the path labeled a, and * 
is the convolution operator. We have chosen to enter 
the relative time delays of the signal explicitly in the 
second term in the sum. In the frequency domain, 
the received spectrum I r (f) takes the form 

I r (f) = J2 G M) W) e~ i2nTaf 

a 

= U/)]TG Q (/) e- a —f 

= Io(f)Gtot(f) (7) 



where all the ugly details of the multipathing are 
contained in the single complex gain Gtot(f) = 
E G «(/) e~ t27rT *f, and the I Q (f) stands cleanly 
apart. Gtot(f) is likely to have strong frequency de- 
pendence, caused by differential phase winding from 
the relative delays. Clearly, application of algorithms 
such as Eqns 5 and 6 will require spectral resolution 
for the frequency channels of Af < l/r max , where 
T m ax is the greatest relative delay. The equivalent 
statement for algorithms that operate in the time do- 
main is that the number of lags used in constructing 
corrective filters must be sufficiently great to capture 
the full range of delays. 
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6. Eigen Decomposition and Gain 
Calibration 

Interferometer arrays measure correlations be- 
tween the signals arriving at the array elements. 
Correlation matrices, such as shown in Fig. 5, can 
be used to summarize the covariances of the signals 
throughout each integration time. Since astronomers 
seek to detect sources that are too weak to detect in 
the short time of a correlator dump interval, they ex- 
tract information from the measured covariances (or 
fringe visibilities) by searching for subtle correlations 
through inversion of the fourier integral of Eqn. 1. 
This is an effective way of reorganizing the correla- 
tion information by restricting the range of interest 
to only those signals that could have originated from 
within the field of view of their telescopes' primary 
beams pb{x,u). Since the Earth's rotation causes 
the u,v values to change over time, the coefficients of 
the covariance matrix are expected to change with 
time; making maps is an efficient way to recover the 
correlations in a physically meaningful way. 

There are other powerful mathematical tools for 
recovering information from correlation matrices, al- 
though these are geared for studying correlations 
over the short integration times for which the coef- 
ficients are expected to remain constant. For exam- 
ple, Leshem, van der Veen and Boonstra (Leshem et 
al 2000 [@] , plus extensive references in Boonstra et 
al 2004 [@]) have extended the eigen decomposition 
methods of signal processing to applications in radio 
astronomy to capitalize on these tools in the area of 
rfi mitigation. Here, we convey some of the more 
intuitive conclusions to the astronomy audience. 

6.1. The covariance matrices 

To begin, Fig. 5 displays a full set of cross power 
spectra for a short observation with the Australia 
Telescope Compact Array - with the spectra po- 
sitioned according to their positions in a correla- 
tion matrix. This data set is "srtcaOl" from Bell 
et al.'s baseband recordings (Bell et al 2001 [@]) 
taken specifically for the purpose of testing rfi miti- 
gation algorithms. The recordings include the time 
sequenced sampling of the outputs from 5 antennas 
with two polarizations, plus two outputs from a dual 
polarized reference antenna aimed at the source of in- 
terference, for a total of 12 input data streams. The 
recorders captured 4 MHz bands with 4 bit preci- 
sion. Correlation, to obtain 512 spectral channels, 



and integration were performed in software (Kocz 
2004 [@]). The rfi affects the central -150 spectral 
channels. 

The total power and cross power spectra in Fig. 5 
illustrate several points: The spectra formed from 
the reference data streams, 1 and 2, in the upper 
left show that the rfi is nearly 100% linearly polar- 
ized, as it appears with equal strength in the total 
and cross power spectra. Further, the cross power 
spectra for 1-2 show that the noise in the passband 
outside the ~150 contaminated channels integrates 
down to a much lower level than in these spectral 
ranges of the total power spectra for cither 1 1 or 
2-2; indeed, this is the motivation for using the cross 
power spectrum instead of a total power spectrum in 
the denominator of the correction term in Eqns. 5 
and 6. Similarly, the total power spectra along the 
diagonal of the matrix in the lower right panels of 
Fig. 5 (involving inputs 3 through 12) scarcely show 
the presence of the rfi in several cases, but the rfi 
clearly appears above the noise in cross power spec- 
tra involving the same antennas. 

6.2. The eigen decomposition 

The eigen decomposition of a 12x12 covariance 
matrix solves for signals that are in common among 
the 12 input data streams. Since there are actually 
512 covariance matrices (one for each of the 512 fre- 
quency channels), the decomposition must be per- 
formed 512 times. In matrix notation, one 12x12 
covariance matrix C, 

/ Ci i Ci 2 ••• C\ 12 \ 
C2 i C22 •■• C2 12 
C3 1 C22 •■■ C312 



C = 



\ C12 1 C12 2 •■• C12 12 / 

decomposes to the product of three 12x12 matrices: 
C = GPG H , where G H is the Hermitian conjugate 
of G and P is a diagonal matrix containing the eigen- 
values. 

/Pi ... \ 
P 2 ... 
P = P 3 ... 

V ... P12 J 

It is convenient to think of the Pi eigenvalues as pow- 
ers. Standard methods for performing the eigen de- 
composition, such as singular value decomposition 
(SVD), produce an ordered list of eigenvalues start- 
ing with the maximum. Should the system have 
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Figure 5. Examples of Power Spectra (diagonal elements) and Cross- Power Spectra (off diag- 
onals) for a 12 signal cross correlation analysis. Each spectrum has 512 spectral channels with 
amplitudes in db. The upper left 2x2 grouping contains the spectra and cross power spectra for 
the reference channels. The lower right 10x10 group has the spectra for the 5 antennas, with 2 
polarizations each, for the ATCA telescopes. The 2x10 and 10x2 groups to the lower left and 
upper right respectively contain the cross correlation spectra that couple the reference antennas to 
the ATCA signals. 



fewer than 12 signals present, the matrix would 
become singular, and some number of eigenvalues 
would be zero. This is not expected of a radio astro- 
nomical system, where there are not only rfi signals 
and celestial sources but also each input data stream 
brings its own noise signal. 

The matrix G containing the eigenvectors 



f 9ii 912 

92 1 92 2 

93 1 92 2 



91 12 

92 12 

93 12 



\ 



\ 912 1 912 2 ■■■ gi2 12 / 



can be visualized as a matrix of complex voltage 
gains gij, giving the coupling of each signal (whose 
voltage is P^ ) to the covariances. Thus, the 
amount of Pj coupled to Cjk is 9ji9kiPi- The top 
row of G specifies the coupling of the 12 Pj to the 
number 1 input, which is the first reference antenna 
in the example here. 

Figure 6 shows the eigenvalue spectra obtained 
from 512 eigen decompositions that were performed 
one spectral channel at a time for the data in Fig. 5. 
The calculation has succeeded in separating the rfi 
and concentrating its power in the first eigenvalue 
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Figure 6. Eigenvalue Spectra for Pt(f) for i=l,12. Ex- 
cept for P\(f), the spectra are very closely spaced, and 
the 11 parallel "passbands" are drawn in alternating light 
and dark colors. 



Pi(/). The other Pi(f) for i > 1 appear to represent 
normal passband power levels expected from receiv- 
ing systems dominated by internal noise. The one 
exception is a narrow rfi spike in P2H) around chan- 
nel 290 that had not been apparent in the spectra in 
Fig. 5. 



6.3. Special case: Array calibration 

An analogy can be made to a common procedure 
in the amplitude and phase calibration of synthesis 
arrays. Typically, an observation is performed on 
a bright, unresolved radio continuum source, whose 
flux density S v is high enough that the source's power 
in the receiving system will dominate over all other 
powers after a short integration time. In the nota- 
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Figure 8. Total power spectra for reference channels 
(1-2) after partial nulling the dominant eigenvalue. 
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Figure 7. Total power spectra for astronomy channels 
(3-12) after partial nulling the dominant eigenvalue. 



Figure 9. Eigenvector phases for the dominant eigen- 
value for all 12 data channels. 
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tion above, this condition means that there is one 
dominant eigenvalue Pi = S v . Under this condition, 
the P matrix degenerates to lxl, and the eigenvec- 
tor matrix G becomes an 12x1 dimensional vector of 
complex gains. In fact, the standard methods of ar- 
ray calibration solve for the "antenna based complex 
gain factors" (for amplitude and phase). These gains 
are a subset of the full eigenvector array, but they en- 
able the reconstruction of the covariance matrix for 
the case with a single powerful source. 
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6.4. Nulling 

Having identified the strength and spectral behav- 
ior of the rfi through its dominance in the eigen de- 
composition of Fig. 6, one is tempted to null the 
rfi by replacing Pi by zero. However, this would 
be slightly too severe a partitioning of the eigenvalue 
matrix, since this would also remove the system noise 
contribution and perhaps distort a weak celestial sig- 
nal contributing noise power to the system. In these 
illustrations reported here, we replace Pi(f) by the 
average of the other Pi(f) (i—2,12), and then recon- 
struct the covariances C(/) from G and the modified 
P matrix. This produces the total power spectra for 
the ten telescope data streams as shown in Fig. 7 and 
the two reference antenna data streams in Fig. 8. 

The partial nulling of the dominant eigen value 
has been effective in recovering total power spectra 
for the telescope inputs (3-12) in which no remnant 
of the rfi is visible. The partial nulling has also re- 
moved the rfi from the reference inputs (1 and 2), 
allowing the second, narrow interferer from P 2 to be 
seen in the spectra for both. Apparently, this narrow 
interferer was only picked up by the reference anten- 
nas and not by the ATCA telscopes themselves. 

The process adopted here to partially null the 
eigenvalue and then recompute an edited covariance 
matrix has the identical effect to the procedure in 
which the unwanted portion of the eigenvalue P ma- 



trix is "projected" back into the covariance matrix 
and then subtracted (e.g. Lcshcm ct al 2000 [@]). 

6.5. Directional information from the eigen 
decomposition 

In the worked example for 12 input data streams, 
the eigen decomposition has identified two interfer- 
ing sources based on the mathematical interpretation 
of the convariance matrix, and it has also enabled 
the nulling of the rfi power. This is accomplished 
without knowledge of location of the interferer, aside 
from having aimed the reference antennas; in fact, 
basing the analysis on the 10 inputs from the ATCA 
dishes alone (and ignoring the reference antenna in- 
puts) leads to a similar result, albeit with a lower 
interference-to-noise ratio in the eigenvalue spectra 
(Kocz 2004 [@]). 

There is directional information contained in the 
phases of the complex eigenvector coefficients. To il- 
lustrate this, Fig. 9 plots the phase spectra for the 
gij coefficients that couple Pi (the dominant eigen- 
value containing the strong rfi power) to each of the 
array antennas. Outside the central channels, the 
phases become random, implying that there is no 
strong broadband signal in these spectral ranges to 
influence these spectral channels. In the center of 
the band where the rfi is strong, the phases show the 
characteristic drift with frequency that is expected 
for interferometer baselines of increasing length (and 
therefore increasing delay) as indicated in Eqn. 7. 
The baseline length to the ATCA telescope number 6 
was much longer that the other baselines between the 
telescopes during this experiment, so the phases for 
data streams 11 and 12 that were produced by tele- 
scope 6 wind very rapidly as a function of frequency 
across the band and therefore are not detected by 
eye in Fig. 9. Radio astronomers might expect to 
decode this phase information to infer directional in- 
formation about the origin of rfi signals - a process 
that would be straight forward if the positions of the 
elements of the array were well known and the el- 
ement gains were isotropic. In practice, the use of 
directional information is hampered by the unknown 
and variable gains of the sidelobes through which the 
RFI enters the antenna (Leshem et al 2000 [@]). 

6.6. Subtleties and cautions 

The eigen decomposition approach shows promise 
for some applications in rfi mitigation. There are 
some subtleties to its application that are covered in 
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greater depth in the series of papers by Leshem ct al 
(e.g., Leshem et al 2000 [@]). (1) Among these is the 
need to prewhiten the noise among the input chan- 
nels in order for the algorithm to perform well; this is 
a step that would require some care in calibration of 
an array. (2) Further, there is the question touched 
upon in section 6.4 where setting the dominant eigen- 
value to zero was attempted in order to null the rfi; 
there is some subtlety in how to choose the right 
amount to remove and how much to leave behind, 
given that the signals from celestial sources are likely 
to appear weakly in the eigenvalues. (3) Astronom- 
ical covariance matrices generally carry signals from 
more sources than the dimensionality of the matrix. 
While rfi might produce the strongest signals, caus- 
ing them to stand out in the eigen decomposition, 
there is also noise signal associated with each of the 
telescope inputs, as well as with the signals from nu- 
merous weak celestial radio sources. Once the num- 
ber of signals exceeds the number of inputs to the 
correlator, a complete recovery of the separate sig- 
nals is no longer possible; in effect, the number of 
unknowns exceeds the number equations available for 
solution. Even the eigen values are then mixtures of 
signals, and the process of nulling runs the risk of 
damaging the astronomical signals. 

7. Summary 

This paper gave an overview of a couple of meth- 
ods for rfi cancellation that show promise for appli- 
cation in radio astronomy observations. There re- 
mains considerable challenge to achieving the dy- 
namic range in rfi cancellation likely to be required. 
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